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Resumen 


El Acuífero Región Carbonífera es uno de los 28 sistemas hidrogeológicos 
existentes en el estado de Coahuila, y también de los que más presiones 
sufre debido a la actividad minera y la sobreexplotación. Este estudio se 
enmarca en la porción oeste de dicho acuífero, abarcando el Ejido Morelos 
y la Colonia El Nacimiento, ambos en el municipio Melchor Múzquiz. En 
dicha área se tomaron muestras en 14 sitios (4 pozos, 6 norias, 2 
manantiales y 2 secciones del río Sabinas), con el objetivo de caracterizar 
desde el punto de vista hidrogeoquímico e isotópico estos cuerpos hídricos 
y así obtener un modelo conceptual de dicho sector. Los resultados 
obtenidos muestran la variabilidad química de las aguas muestreadas, 
clasificándose como aguas de tipo bicarbonatada-cálcica, bicarbonatada- 
sódica, sulfatada-cálcica y clorurada-sódica, lo cual sugiere la ocurrencia 
de variados procesos que afectan la calidad del recurso, entre ellos la 
salinidad. Desde el punto de vista isotópico se propone un origen 
meteórico para la mayoría de los cuerpos hídricos en estudio, cuya zona 
de recarga proviene de la sierra de Santa Rosa. Sin embargo, procesos 
como la evaporación superficial, la disolución de sales, y el intercambio 
iónico generan variabilidad química y fraccionamiento isotópico, lo que 
condiciona firmas isotópicas más pesadas, y aguas más enriquecidas y 
salinas. Estos elementos se representan en un modelo geohidrológico e 
isotópico conceptual, que integra ambas interpretaciones y muestra la 
complejidad hidrogeológica de la región en términos de calidad del 


recurso. 


Palabras clave: hidroquímica, isótopos estables, acuífero Región 


Carbonífera. 
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Abstract 


The Region Carbonifera Aquifer is one of the 28 existing hydrogeological 
systems in the State of Coahuila and one of those that suffer the most 
pressure due to mining activity and overexploitation. This study is framed 
in the western portion of the aquifer, encompassing the Ejido Morelos and 
the Colonia El Nacimiento, both in the Melchor Muzquiz municipality. In 
this area, 14 sites were sampled (4 wells, 6 norias, 2 springs, and 2 
sections of the Sabinas River) with the aim to characterize these water 
bodies from the hydrogeochemical and isotopic point of view, thus 
obtaining a conceptual model of this sector. The results obtained show 
the chemical variability of the sampled waters, classifying waters as 
bicarbonate-calcium, bicarbonate-sodium, sulfated-calcium, and chloride- 
sodium, which supports the occurrence of various processes that affect 
the quality of the resource, among them the salinity. From the isotopic 
point of view, a meteoric origin is proposed for most water bodies under 
study, whose recharge zone comes from Santa Rosa mountains. However, 
surface evaporation, the dissolution of salts, and ¡onic exchange 
processes generate chemical variability and isotopic fractionation, which 
conditions heavier isotopic signatures and more enriched and saline 
waters. These elements are represented in a conceptual hydrochemical 
and isotopic model, which integrates both interpretations and shows the 
hydrogeological complexity of the region in terms of the quality of the 
resource. 


Keywords: Hydrochemistry, stable isotopes, Región Carbonifera Aquifer. 
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Introducción 


La salinidad en acuíferos mexicanos es un tema preocupante en las 
últimas décadas. De acuerdo con Raynal-Gutiérrez (2020), en la república 
mexicana se han identificado 32 acuíferos con problemas de salinidad, 
principalmente en climas áridos y semiáridos, donde prevalecen las bajas 
precipitaciones y altas tasas de evaporación. Además, en al menos 18 
acuíferos costeros se reportan problemas de intrusión marina, la mayoría 
en la península de Baja California y en Sonora (Torres-Martínez, 
Mahlknecht, Hernández-Antonio, 8 Mora, 2017; Raynal-Gutiérrez, 2020). 


En otras partes del país la salinidad se ha incrementado por niveles 
de aguas subterráneas poco profundas, disolución de minerales 
evaporíticos, la presencia de aguas connatas O fósiles, y la 
sobreexplotación de los recursos hídricos (Wurl, Gamez, Ivanova, 
Lamadrid, 8 Hernandez-Morales, 2018; Gilabert-Alarcón et al., 2018). 
Sumado a esto, los efectos del cambio climático indican problemas futuros 
en la disponibilidad y calidad del recurso hídrico, como así lo refieren 


Hassani, Azapagic y Shokri (2021) en su modelo predictivo global. 


En el estado de Coahuila, los recursos hidrogeológicos están 
distribuidos en 28 acuíferos que abarcan diferentes climas, condiciones 


geológico-estructurales y diferentes presiones de explotación. De hecho, 
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a nivel de estado, datos compilados en 2018 muestran que un no 
despreciable 32 % de los acuíferos sufre fenómenos de salinidad mientras 
que siete acuíferos presentan sobreexplotación (Conagua, 2020). El 
Acuífero Región Carbonífera (ARC) no está exento a los efectos de 
salinización debido a una compleja geología y a la predominancia de un 
clima árido/semiárido extremo. En todo el acuífero son limitadas las 
investigaciones relacionadas con la composición química y la calidad del 
agua de los cuerpos hídricos que lo componen. La Comisión Nacional del 
Agua (Conagua) es la que con algún grado de detalle ha estudiado el 
acuífero a partir de su informe “Actualización de la Disponibilidad de Agua 
en el Acuífero Región Carbonífera, estado de Coahuila” (Conagua, 2020), 
donde detalla los aspectos más relevantes de este sistema en cuanto a la 
geología, geomorfología, características hidrogeológicas, y los cálculos de 
la disponibilidad y balance hídrico. Un estudio más detallado lo presentan 
Batista-Cruz et al. (2022), quienes estudiaron 155 muestras de aguas 
tomadas en pozos, norias, manantiales y fuentes superficiales. El estudio 
pormenoriza aspectos importantes de la calidad del agua con base en la 
NOM-127-SSA1-1994 (DOF, 2000) para fines de consumo y, a Su vez, 
muestra la compleja situación del acuífero, cuya calidad es cuestionable 
para actividades agrícolas. Sin embargo, debido a la variabilidad química 
observada por dichos autores, los estudios hidrogeoquímicos e isotópicos 
de áreas pequeñas dentro del acuífero pueden ayudar a comprender 
mejor qué procesos deterioran la calidad del recurso hídrico y, en tal 
sentido, cómo enfrentar en municipios y ejidos la problemática del abasto 


de agua para diversos fines. 
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Por su parte, los estudios isotópicos (p. ej., isótopos estables de 
oxígeno y deuterio) son prácticamente nulos, con solo información en los 
acuíferos cercanos de Cuatrociénegas y el Hundido (Johannesson, 
Alejandra, €: Kilroy, 2004; Felstead, Leng, Metcalfe, 8 Gonzalez, 2015). 
La sinergia entre química e isotopía permite identificar disímiles, y 
complejos procesos hidrogeológicos y ambientales, y entender aquellas 
causales que pueden modificar la composición química prístina de un 


determinado cuerpo hídrico, así como su origen. 


El presente estudio tiene como objetivo caracterizar, desde el punto 
de vista hidrogeoquímico e ¡isotópico, las aguas superficiales y 
subterráneas del Ejido Morelos y la Colonia El Nacimiento en el municipio 
de Melchor Múzquiz, Coahuila, para identificar los procesos geo- 
hidrogeológicos e isotópicos que sufren las aguas subterráneas en su 
movimiento por el medio natural y, a su vez, elaborar con ellos un modelo 
conceptual del origen de las aguas en este sector del Acuífero Región 


Carbonífera. 


Área de estudio 


El área de estudio se centra en el Ejido Morelos y la Colonia El Nacimiento, 
dos zonas contiguas que cubren un área total de 521 km?, localizadas en 
el municipio Melchor Múzquiz y en la porción oeste del ARC (Figura 1). El 
acuífero, clasificado como de tipo libre, está constituido en su parte 
superior por sedimentos aluviales y conglomerados polimícticos de la 
Formación Sabinas. En su parte intermedia se presenta la Formación 


Olmos y a mayor profundidad las rocas calcáreas —fundamentalmente de 
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las Formaciones Aurora y Cupido del Cretácico Inferior— constituyen la 
unidad acuífera que puede presentar condiciones de confinamiento debido 
a que está sobreyacida por rocas arcillosas como lutitas y limolitas (SGM, 
2003; Conagua, 2020). 
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Figura 1. Mapa de ubicación geográfica del área de estudio dentro de 
los límites del Acuífero Región Carbonífera en el estado de Coahuila. La 
base geológica fue tomada de SGM (2003). En el mapa: M-municipio 


Melchor Múzquiz. 
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Dentro de los límites del acuífero, la profundidad del nivel estático 
varía de 5 a 30 metros. Sin embargo, dicha profundidad puede alcanzar 
hasta 50 metros en algunos pozos perforados en calizas que se localizan 


en la sierra de Santa Rosa. 


La hidrología superficial se caracteriza por cauces intermitentes y 
efímeros, teniendo muy pocos cauces permanentes, como el caso del río 
Sabinas, que transita por el área de estudio y cuyos datos isotópicos se 


presentan en este estudio. 


Materiales y métodos 


Muestreo y técnicas analíticas 


Los trabajos de muestreo se llevaron a cabo en 14 sitios (4 pozos, 6 
norias, 2 manantiales y 2 secciones del río Sabinas) (Tabla 1). En cada 
sitio se tomaron las coordenadas y la altitud, y se midieron los parámetros 
físico-químicos (pH, conductividad y sólidos disueltos totales) (Tabla 1) 
con la ayuda de un equipo multiparamétrico marca HANNA HI98194., Al 
terminar cada medición, el equipo se lavó con agua desionizada para 
evitar errores de medición, así como para su propio cuidado y buen uso. 
Para el caso de los muestreos químicos de elementos mayoritarios (Cl, 
SO4?", HCO>", Na*, Ca?*, K*, Mg?*) se colectaron 500 mililitros de muestra 
de agua en botellas de polietileno de alta densidad (HDPE), previamente 


lavadas en una solución de ácido nítrico. Además, las muestras de agua 
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para cationes se acidificaron en campo con HNO3a pH = 2. Para el caso 


de los análisis isotópicos (9180 y 9D), las muestras de agua fueron 


extraídas con una jeringa de 50 ml acoplada a un filtro VWR cónico estéril 


de 25 mm de diámetro y 0.45 micrómetros de tamaño de poro. Todas las 


muestras se tomaron por duplicado y al finalizar cada muestreo fueron 


refrigeradas en hieleras herméticas a una temperatura entre 1 y 5 grados 


Celsius para evitar 


cualquier tipo de contaminación. 


Todos los 


procedimientos de muestreo se basaron en los trabajos de Nollet (2000). 


Tabla 1. Datos físico-químicos de las muestras de agua en estudio. 


usd Tipo 5H SDT C.E. HCOz" | SO4?- (e Na* k+ Ca?+ Mg?* |Error 
(ppm) |(1S/cm) | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (%) 
Casa Blanca Pozo 8.08 | 1579 3 160 538 24.2 285.4 307 1.31 22.3 2.4 8.34 
129 Pozo 8.3 385 460 237.9 21.5 12.3 3.7 0.4 81 12.2 |-5.28 
124 Pozo 7.8 1319 4 770 2 446.1 20.6 397.6 1244 5.1 37 18.5 |-5.38 
Pozo La Flor Pozo 7.87 855 1711 476 71.2 26 188 1.23 12.3 2.5 5.19 
143 Noria 7.7 451 540 305 11 10.01 3.7 0.4 101 8.8 -3.72 
135 Noria 7.6 565 650 353.8 31.2 14.2 6.2 1.6 140 3.6 -5.16 
13 Noria 7.9 176 320 109.8 11 10.7 4,4 1.6 32 6.4 -0.54 
17 Noria 7.6 472 730 79.3 194,3 60.4 46 3.9 74 14.4 0.51 
Ejido Morelos Noria 6.69 895 1790 345 38.2 81.9 57.9 1.58 127 10.1 |-5.23 
Iglesia Nacimiento Noria 7.52| 2146 4 305 499 9.6 474.4 383 1.8 20.6 3.4 9.42 
Parquecito Río 6.79 551 1103 351 42.2 7.8 4.7 0.73 117 7.9 1.03 
Kikapoo Río 8.03 | 383 766 185 53.7 9.7 10.1 0.99 69 10.8 |-4.04 
Ojito Manantial | 7.33 589 1177 409 24.9 11.6 9.71 2.63 118 0.96 7.79 
137 Manantial | 7.8 357 450 231.8 20.6 14.2 2.3 0.4 81 7 -1.07 
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Los análisis físico-químicos a escala de laboratorio se desarrollaron 
en el Laboratorio de la Facultad de Química de la Universidad Autónoma 
de Coahuila (UADEC), mediante diversos métodos volumétricos y 
analíticos respaldados por las normas mexicanas correspondientes. Para 
el caso de los cationes Ca*+ y Mg?*, sus concentraciones se calcularon 
multiplicando la dureza temporal, calculada previamente por titulación 
con EDTA 0.2N (SCT, 2001c) por el factor 0.4 y 0.24, respectivamente. 
Por su parte, las concentraciones de Na* y K* se obtuvieron mediante 
espectrofotometría de absorción atómica con llama (AA240FS Varian) 
utilizando curvas de calibración con cinco estándares de 0.4 a 1.2 mg/l y 
de 0.3 a 1.5 mg/l, respectivamente, a partir de estándares certificados de 
1 000 mg/l (SCFI, 2001b). Las concentraciones de los aniones HCO3x', 
SOa4?- y Cl" se adquirieron por titulación con H2SO4 a 0.02N (SCFI, 2001a), 
BaClz a 0.1N (SCFI, 2014), y AgNO3 a 0.01N (SCFI, 2001d), 


respectivamente. 


Los análisis isotópicos se hicieron con el analizador isotópico de 
agua Picarro L2130-1, propiedad del Centro de Investigación Científica y 
de Educación Superior de Ensenada Baja California. El equipo, basado en 
la tecnología “Cavity Ring-Down Spectroscopy” (CRDS) mide la razón 
isotópica de 9D, 9180 con alta precisión (errores bajos inferiores a 0.1 %o) 
(Picarro Inc., 2008). Las muestras filtradas y sin ácidos adicionados son 
contenidas en pequeños viales (hasta 200 pul) y colocadas en un 
automuestreador del equipo, conjuntamente con agua desionizada y los 
respectivos estándares internos Pi-Z (9180 = 0.3 %o; 9D = 1.8 %o); Pi-M 
(9180 = -20.6 %o; 9D = -159.0 %o) y Pi-D (9180 = -29.6 %o; 0D = - 
235.0 %o); y los estándares internacionales certificados VS9MOW2 (9180 
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= 0 %o; 0D = 0 %o), GISP (9180 = -24.7 %o; 0D = -189.9 %o) y SLAP 
(9180 = -55.5 %o; 0D = -428.0 %o). En cada análisis se obtienen 10 
mediciones, donde las primeras cinco son ignoradas (para eliminar el 
efecto memoria), y las últimas cinco se consideran para la corrección de 
los datos sin procesar. La corrección de los datos se realiza con base en 
los estándares internacionales e internos por el método de regresión 


lineal. 


La reproducibilidad de los análisis para 9180 y 9D fue de alrededor 
de + 0.1 %o y + 0.4 %o, respectivamente. Esos datos son reportados en 
comparación con el estándar Vienna Standard Mean Ocean Water 2 (V- 
SMOW2), mientras que los resultados se informaron en notación 9, como 
desviaciones de las relaciones isotópicas en partes por mil (%o) de 


acuerdo con los criterios de Craig (1961). 


Resultados 


Las 14 muestras en estudio presentan diferencias físico-químicas entre sí 
(Tabla 1), lo que presupone la ocurrencia de variados procesos 
hidroquímicos en dichos sitios. El error de balance calculado a partir de la 
sumatoria de aniones y cationes en meq/l arrojó valores inferiores al 10 


%, lo que indica una calidad óptima de los datos (Appelo 8, Postma, 1994). 


En detalle, con la excepción del pozo 129, todas las muestras 
provenientes de los pozos presentan una mayor concentración en sólidos 
disueltos totales (SDT) y en la conductividad eléctrica (C.E.) (Figura 2A). 
En detalle, los pozos Casa Blanca (1 579 ppm) y 124 (1 319 ppm), y la 
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noria Iglesia Nacimiento (2 146 ppm) tienen concentraciones de SDT por 
encima del valor máximo permisible estipulado por la Norma Oficial 
Mexicana NOM-127-SSA1-1994 (Figura 2A; SDT-1000 ppm), por tanto 
presentan una deficiente calidad de agua para fines de consumo humano. 
El resto de las norias, fuentes superficiales y manantiales muestran 


valores de SDT inferiores a 895 ppm y C.E. por debajo de 1 177 uS/cm. 
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Figura 2. Análisis químico comparativo de las muestras de agua: A) 
contenido de sólidos disueltos totales y conductividad eléctrica; B) 
concentraciones de los elementos mayoritarios en las aguas. C.E.: 


conductividad eléctrica; SDT: sólidos disueltos totales. 


42 
6 = 2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 28-69. DOI: 10.24850/j-tyca-2024-03-02 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


En relación con la composición química de las aguas (Figura 2B), se 
observa que las muestras más enriquecidas en los constituyentes 
mayoritarios son nuevamente los pozos 124 y Casa Blanca, así como la 
noria Iglesia Nacimiento. Los mayores contenidos en Cl se reportan de la 
siguiente manera: Iglesia Nacimiento (2 146 mg/l); Pozo Casa Blanca (1 
579 mg/l), y Pozo 124 (1 319 mg/l), mostrando un enriquecimiento de 
8.5, 6.3 y 5.2 veces, respectivamente, en comparación con el límite 
máximo permisible que establece la Norma Oficial Mexicana NOM-127- 
SSA1-1994 (CI = 250 mg/l). Para el elemento Na* se reportan 
enriquecimientos de hasta seis veces en la muestra del pozo 124 con 
respecto a la norma de calidad de agua antes mencionada. De forma 
general, los cuerpos superficiales y los manantiales presentan 
concentraciones bajas en la mayoría de los constituyentes químicos lo 


cual demuestra que no están siendo afectados por procesos profundos. 


Desde el punto de vista isotópico, las aguas en análisis varían entre 
-3.29 %o y -7.36 %o para el isótopo 9180; mientras que valores entre - 
23.02%o y -48.24 %o se reportan para el deuterio (Tabla 2). De forma 
general se observa que las aguas superficiales (p. ej., Kikapoo y 
Parquecito) presentan los valores más empobrecidos (más ligeros) en 
ambos isótopos en relación con el resto de los cuerpos de agua, mientras 
que la muestra más enriquecidas en ambos isótopos (más pesadas) es la 
Noria (9180 = -3.29 %o). El pozo Aparicio no se encuentra en el área de 
estudio, pero debido a su localización en la sierra de Santa Rosa puede 
ayudar a comprender la relación de esta cadena montañosa con la 


hidrología y los procesos de recarga ocurridos en la región. 
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Tabla 2. Composición isotópica de oxígeno y deuterio de las aguas en 


estudio (msnm: metros sobre el nivel medio del mar). 


Muestras Descripción ys Re a 
(%o) (Y%o) | (msnm) 
Kikapoo Río -7.36 -48.24 513 
Parquecito Río -6.66 -43,49 455 
Ejido Morelos Noria -6.09 -38.83 492 
Iglesia Nacimiento Noria -4.79 3.73 489 
Noria Noria -3.29 -23.09 448 
Pozo La Flor Pozo -6.24 -39.83 513 
Pozo Casa Blanca Pozo ES -34.87 605 
Pozo Aparicio Pozo -7.04 -44.62 605 
Ojito Manantial -5.02 -32.85 564 
Divi Manantial 100 -34.21 570 
Agua meteórica Múzquiz | Agua meteórica 7.23 -47.8 500 


Fuente: La composición del océano fue tomada de Craig (1961). 
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Discusión 


Procesos hidrogeoquímicos e isotópicos que controlan 
la variabilidad química en los cuerpos hídricos 


En la Figura 3 se presenta el mapa hidroquímico regional basado en los 
diferentes diagramas de Stiff (1951), lo que ofrece una visión general de 


los diferentes tipos de agua y su localización espacial. 
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Figura 3. Mapa del quimismo de las aguas del área de estudio a partir 
de los diferentes gráficos de Stiff. En el mapa: cubos negros-pozos, 
triángulos verdes-norias, círculos rojos-manantiales y estrellas azules- 


ríos y arroyos. 
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Como se indica en la Figura 3, se identifican nueve muestras de 
agua de tipo bicarbonatada-cálcica y/o magnésica; tres muestras de tipo 
bicarbonatadas-sódicas; una muestra de tipo clorurada-sódica y/o 
potásica, y una muestra clasificada como sulfatada cálcica y/o magnésica. 
En el caso de la noria Ejido Morelos, el patrón del diagrama de Stiff, 
aunque mayoritariamente de tipo bicarbonatado cálcico, tiene 
incrementos de cloruros y sodio, lo que sugiere procesos de aumento en 
la salinidad. De forma general, a pesar de ser un área no tan extensa, 
con una superficie total de 521 km?, se observa gran variabilidad química 
en los diferentes cuerpos de agua, lo cual es el resultado de variados 
procesos geohidrológicos y ambientales. En tal sentido, el diagrama de 
Piper (1944), modificado por Kelly (2006), ofrece una interpretación más 


completa. 


Kelly (2006) definió los principales procesos hidroquímicos 
responsables de la variabilidad química en los cuerpos de agua, aportando 
un enfoque diferente al rombo de Piper (1944), en el que establece siete 
nuevos campos (Figura 4). En dicha figura se observa que las aguas 
clasificadas como bicarbonatadas-cálcicas (p. ej., muestras 129, 135, 
137, Ojito, 13, Parquecito, 143 y Kikapoo) se localizan en el campo de 
aguas dulces, por lo que desde el punto de vista de calidad de agua son 
óptimas para diversas actividades. La composición de tales aguas 
proviene de la disolución de minerales presentes en rocas carbonatadas 
ricas en Ca?*, Mg?* y HCO3" (Choi et al., 2014), que conforman la sierra 
de Santa Rosa en las formaciones La Peña y Aurora. Sin embargo, se 


observa que la muestra de la noria Ejido Morelos se aleja del vértice, con 
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una tendencia hacia el campo definido por Kelly (2006) como “mezcla 


conservativa ligera” por intrusión salina, como se observa en la Figura 4. 


Simbología Ca/Mg-SO, M Pozos 
Agua dulce e y A Norias 


Intrusión marina ligera yo O Manantiales 
Evolución Natural ligera A / *k Río 
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Figura 4. Diagrama Kelly (2006) que muestra los principales procesos 


hidroquímicos que controlan el quimismo de las aguas. 
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Las aguas del tipo bicarbonatadas-sódicas son reportadas para los 
pozos La Flor, 124 y Casa Blanca (Figura 4). El aumento en la 
concentración de sodio en estas aguas parece estar relacionado con 
procesos de interacción agua-roca donde existe elevado intercambio 
iónico, probablemente entre arcillas, fuente importante de suministro de 
sodio. Dichas secuencias podrían estar en profundidad interactuando con 
los flujos del acuífero. Para el caso del pozo Casa Blanca y la noria Iglesia 
Nacimiento se observa una tendencia hacia el vértice de las aguas 
cloruradas sódicas, lo cual es un indicio de salinización de dichas aguas. 
En este caso y, por tratarse de pozos y norias, la disolución de sales 
existentes en el medio geológico, probablemente de naturaleza 
evaporítica muy comunes en el norte de Coahuila (Padilla-y-Sánchez, 
2007; González-Sánchez, Puente-Solís, González-Partida, 8 Camprubí, 
2007; Eguiluz-de-Antuñano, 2001), podría justificar dicho proceso y, por 
ende, el aumento en las concentraciones de Cl, Na*, SDT y la 
conductividad eléctrica discutidos previamente. No obstante, los procesos 
superficiales, como la evaporación o la infiltración de aguas o soluciones 
muy evaporadas, también podrían inducir sales al medio, ya que se trata 
de zonas semiáridas donde la evaporación es dominante sobre las 


precipitaciones (Batista-Cruz et al., 2022). 


Una sola muestra (p. ej., noria 17) es clasificada como sulfatada- 
cálcica y/o magnésica, y se localiza en el campo de intrusión salina (Kelly, 
2006). Aguas de este tipo pueden deber su composición a diferentes 
procesos, entre los que se destacan los siguientes: a) lavado de terrenos 


formados en ambiente marino; b) oxidación de sulfuros que se 


49 
2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 28-69. DOI: 10.24850/j-tyca-2024-03-02 


(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 1) Check for updates 
OPEN ACCESS 


Tecnología y 

CienciaszAgua 
encuentran ampliamente distribuidos en rocas ígneas y sedimentarias, y 
c) descomposición de sustancias orgánicas. En el área de estudio, los 
procesos a) y b) son factibles debido a que, como ya se expresó, la 
geología del estado presenta secuencias evaporíticas de las antiguas 
condiciones marinas restringidas que imperaron en la cuenca Sabinas 
desde el Triásico hasta el Cretácico Superior (Eguiluz-de-Antuñano, 
2001), a la cual pertenece el área de estudio. La existencia de secuencias 
geológicas con rocas y minerales solubles como el yeso, anhidrita, halita, 
y sales magnésicas y potásicas es ampliamente documentado en toda la 
región, lo cual representa un problema para la calidad del agua en varios 
municipios por salinización (Batista-Cruz et al., 2022). Sin embargo, la 
oxidación de sulfuros también es una vía no descartable de aporte de 
sulfato a las aguas. En este sentido, la sierra de Santa Rosa y gran parte 
de la geología coahuilense posee amplia distribución de depósitos 
minerales metálicos y submetálicos con importantes minas galena y 
hierro en explotación, por lo que ambos procesos parecen ser factibles 


como aportadores de sulfato a las aguas. 


A partir del diagrama de Mifflin (1988) (Figura 5) se identificaron 
los diferentes tipos de flujos subterráneos con base en la concentración 
de los iones Na*, K*, Cl" y SO4?”. Como se observa en el gráfico, las 
muestras en estudio se localizan en los tres flujos identificados por el 


autor (local, intermedio y regional). 
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Figura 5. Diagrama binario de Mifflin (1988), indicando la procedencia 
de los flujos local, intermedio y regional. En la figura se detallan además 
los diferentes tipos de aguas identificadas en los diagramas de Stiff y 


Piper para una mejor comprensión. 


Las aguas pertenecientes al flujo local (manantiales 137 y Ojito; 
pozo 129; norias 13, 143, 135; y las muestras de aguas superficiales 
Kikapoo y Parquecito) corresponden a aguas poco evolucionadas y de 
reciente infiltración, que coinciden con una naturaleza química 


bicarbonatada-cálcica/magnésica y con los valores más bajos de las 
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relaciones Na++K* y CI+SO4*. La muestra de Ejido Morelos, que 
anteriormente se habían identificado en el diagrama de Stiff y Piper con 
algo de salinidad debido a la intrusión salina, corrobora su tendencia hacia 
aguas más mineralizadas que se localizan en el flujo intermedio (Figura 
5). En este campo se localizan además las muestras de agua provenientes 
de los pozos La Flor y Casa Blanca, así como las norias 17 e Iglesia 
Nacimiento. Estas aguas presentan una mayor evolución 
hidrogeoquímica, dada por un mayor tiempo de residencia en el medio 
geológico y la profundidad a la cual se infiltraron. De hecho, las aguas de 
Iglesia Nacimiento van evolucionando hacia la presencia de flujos 
regionales. De forma general, todos estos cuerpos hídricos tienen 
mayores concentraciones químicas, como lo demuestran los elevados 
valores en Na++K* y CI+SO4”, en comparación con el resto de las 
muestras (Figura 5). De acuerdo con Chapelle (2000), los sistemas de 
flujo intermedio corresponden por lo general a acuíferos confinados de 
profundidad moderada (< 300 metros) y que se diferencian de los 
sistemas de flujo locales sobre todo en la extensión de la conexión con la 
superficie y con la zona no saturada. Al tratarse de pozos y norias, la 


presencia de estos flujos debajo de estos sitios parece estar justificada. 


Por último, las aguas del pozo 124, que son las más enriquecidas 
en Nat y HCOz" (ver también Tabla 1 y Figura 4), se localizan 
químicamente dentro del flujo regional, por lo que corresponden a una 
zona de depositación de solutos, donde el agua se va saturando de iones, 
producto a la mineralización a lo largo de la trayectoria recorrida por el 
subsuelo. Así, se consideran las aguas más viejas, maduras y 
evolucionadas (Mifflin, 1988). 
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Composición isotópica de las aguas 


En estudios hidroquímicos, el uso de los isótopos estables (átomos de un 
mismo elemento con igual número atómico y diferente número másico) 
han influido de forma positiva en la interpretación de los procesos que 
forman parte del ciclo hidrológico y, por ende, permiten clarificar el origen 
del agua y los procesos de fraccionamiento que afectan su composición 
(Clark 8 Fritz, 1997; Mook, 2001). Con el objetivo de integrar la 
información isotópica con la composición hidrogeoquímica, en esta 
sección se presenta el origen de las aguas y el modelo conceptual del área 
de estudio. Lamentablemente no se cuenta con datos isotópicos de todas 


las muestras, por lo que el enfoque debe tomarse como general. 


El gráfico isotópico oxígeno versus deuterio de la Figura 6 muestra 
los diferentes procesos de fraccionamiento presentes en los cuerpos 
hídricos en estudio. Las muestras superficiales de Kikapoo (0180 -7.36 
%o; 0D -48.24 %o) y Parquecito (9180 -6.66 %o; 9D -43.49 %o); los 
pozos Aparicio (9180 -7.04 %o; 9D -44.62 %o) y La Flor (9180 -6.24 %o; 
9D -39.83 %o), así como la noria Ejido Morelos (9180 -6.09 %o; dD - 
38.83 %o) tienen un origen meteórico, procedente de las precipitaciones 
locales descargadas en la sierra de Santa Rosa (9180 -7.23 %o; 90D -47.8 
%o), como lo muestra su localización sobre la línea mundial de aguas 
meteóricas de Craig (1961). Los ligeros cambios isotópicos hacia valores 
más pesados observados en el Pozo La Flor y la noria Ejido Morelos en 
relación con el agua meteórica se deben probablemente a diferencias del 
altitud y a procesos de interacción agua-roca o intercambio iónico 
(Giggenbach, 1991). 
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Figura 6. A) Gráfico isotópico oxígeno versus deuterio donde se han 
representado la línea mundial de aguas meteóricas y la composición 
isotópica del océano (VSMOW), ambas tomadas de Craig (1961); B) 


gráfico 9180 versus cloruros. 
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Los cuerpos hídricos con las firmas isotópicas más pesadas que las 
aguas meteóricas locales son Divi (3180 -5.38 %o; 90D -34.21 %o); Casa 
Blanca (9190 -5.33 %o; 0D -34.87 %o); Ojito (980 -5.02 %o; 9D -32.85 
%o); Iglesia Nacimiento (9180 -4.79 %o; 9D -31.73 %o), y Noria (9180 -— 
3.29 %o; 90D -23.09 %o). En un principio, estas aguas deben provenir de 
las precipitaciones locales, pero han sido fraccionadas por procesos de 
evaporación superficial y/o disolución de minerales existentes en el 
subsuelo. Debido a la intensa evaporación superficial y las escasas 
precipitaciones de la región, al igual que la limitada profundidad de las 
estas obras (no superan los 50 metros), estos procesos parecen ser 
factibles. La Figura 7 y Figura 8 muestran las particularidades geológicas, 
hidrológicas y ambientales que parecen controlar el quimismo de las 
aguas en el área de estudio. Desde el punto de vista geológico afloran 12 
formaciones con edades que van del desde el Cretácico Inferior hasta el 
Neógeno y dos litologías cuaternarias innominadas del Holoceno 
compuestas por depósitos sedimentarios aluviales (SGM, 2003). El origen 
de la salinidad y el enriquecimiento en algunos ¡ones como el sulfato 
podría estar relacionado con las condiciones de depósito y los ambientes 
sedimentarios que dieron lugar a dichas formaciones geológicas. Como se 
muestra en la Figura 7, en la sierra de Santa Rosa, el anticlinal Santa Ana 
y anticlinal Loma Gorda, las formaciones son de origen sedimentario, con 
predominancia de calizas y lutitas depositadas en ambientes marinos 
desde profundos (formación La Peña) hasta someros (formaciones 
Cupido, Aurora) y costeros de tipo estuarino tipo sabkha (formación La 


Mula) (Figura 7). En la base de las sierras y las mesetas afloran las 
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formaciones Buda, Eagle Ford, Austin, Upson, San Miguel y Olmos, cuya 
composición petrológica es mayormente de tipo lutitíco, limolítico y 
arcilloso, con presencia también de calizas. Los ambientes de 
sedimentación varían desde neríticos hasta transicionales de tipo deltaico 


o prodelta. 
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Sierra de Santa Rosa 


Santa Ana 


Ichor Muzquiz 


E 


Ambientes sedimentarios 
Holoceno [_Qhoal | Depósitos aluviales Depósitos aluviales (continental aluvial) 
Holoceno | QpthoCgp| Conglomerado poligenético Conglomerado poligenético (continental aluvial) 
Neógeno [ TpiCgo ] Fm Sabinas Fm Sabinas (continental fluvial) 
Fm Olmos Fm Olmos (marino somero y deltaico) 


[ KcAr-Lu | Fm San Miguel Fm San Miguel (deltaico) 
| KcLu-Lm | Fm Upson Fm Upson (plataforma de aguas someras y prodelta) 


Fm Austin (nerítico con aporte terrígeno) 
KcetLu-Cz| Fm Eagle Ford Fm Eagle Ford (marino somero) 


| KceCz | Fm Buda Fm Buda (nerítico 
m Del Río | [Kacecztu] Grupo Westilta Fm Del sl a 
Fm Aurora Fm Aurora (plataforma somera) 
Fm La Peña Fm La Peña (aguas profundas en zonas batiales) 
| KapCz | Fm Cupido Fm Cupido (plataforma somera) 


| KhbLu-Cz | Fm La Mula Fm La Mula (marinas a estuarinas tipo sabkha) 


Cretácico superior 


Cretácico inferior 


Figura 7. Mapa geológico del área de estudio realizado con la base 
geológica 1:250 000 del SGM (2003). En la leyenda del mapa se 
presentan las formaciones existentes en el área y se sintetizan los 
ambientes de sedimentación de acuerdo con el léxico estratigráfico 
mexicano tomado de SGM (2003). La litología depósitos aluviales 


(Qhoal) está representada sin color de fondo en el mapa. 
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A Nx > 
Anticlinal | 
Santa Ana 


Evaporación PAP LF Río Sabinas EM, IN, CB N,D,O 


" (6180 -7.04 %o) ., (80 -6.24 %o) ., (60 -7.36 a -6.66 %o) , (80 -6.09 a -4.79 %o) (680 -5.02 a -3.29 %o) 
7 Escurrimientos (6D -44.62%o) (6D -39.83%o) (5D -48.24 a -43.49%o0) (8D -38.83 a -31.73%o0) ” (6D -34.87 a -23.09%o) 


A superficiales 


Figura 8. Modelo geohidrológico e isotópico conceptual del área de 
estudio representando los principales procesos hidrológicos y 
ambientales que controlan el quimismo de las aguas. Las letras en la 


leyenda inferior están referidas en el texto. 


Numerosos trabajos se han referido a la presencia de sales y 
minerales, como el yeso, en algunas de estas formaciones en varias 


partes de Coahuila. De acuerdo con Santamaría et al. (1991), en la 
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formación La Mula, por ejemplo, existen horizontes de minerales solubles 
como anhidrita y yeso. Esta formación aflora en nuestra área de estudio, 
lo cual sugiere que los flujos superficiales y subterráneos de la sierra de 
Santa Rosa podrían lixiviar dichos minerales, y aportar Ca?* y SO” a las 
aguas, como sucede en la muestra noria 17. Este proceso ha sido 
reportado previamente en algunas sierras que rodean el área de 
Cuatrociénegas (Mamer 8: Newton, 2017) y en la Sierra Menchaca (Czaja, 
Estrada-Rodríguez, 8: Flores-Olvera, 2014). 


Para el caso de las evaporitas, la existencia de formaciones 
formadas en ambientes deltaicos, estuarinos y marinos someros 
semirrestringidos (Figura 7) podrían contener minerales como la halita 
(NaCl) en sus estratos o secuencias. Sin embargo, aunque en teoría 
parece ser factible, al menos en la sierras que rodean el área de estudio 
no hay reportes previos u observaciones de campo sobre la presencia de 
este mineral u otros minerales aportadores de Cl: y Na* a las aguas, por 
lo que parece tratarse de procesos superficiales. De hecho, tampoco se 
reportan sales en los valles aluviales formados por la formación Sabinas, 


los conglomerados poligenéticos y los depósitos aluviales. 


Los procesos superficiales de carácter ambiental se presentan en el 


modelo de la Figura 8. 


El modelo comienza con la caída de aguas meteóricas sobre las 
sierras que rodean la subcuenca Esperanzas, donde se localiza el área de 
estudio (Figura 8). Los procesos de recarga-descarga permiten que los 
valles y planicies obtengan volúmenes importantes de agua y con ello las 
firmas isotópicas iniciales de tipo meteórica (9180 -7.23%o; 9D -47.8%o) 


se van modificando hacia valores ligeramente más pesados en ambos 
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isótopos, como se evidencia en las muestras de los pozos PAP (09180 - 
7.04%o; 9D -44.62%o) y LF (9180 -6.24%o; 9D -39.83%o). Estas firmas 
son consistentes con procesos de infiltración en las secuencias 
carbonatadas e intercambio iónico. Por su parte, la composición isotópica 
del río Sabinas (9180 -7.36 a -6.66%o; 0D -48.24 a -43.49%o), tomada 
de las muestras Parquecito y Kikapoo, está en el rango de la composición 
del agua meteórica, lo cual confirma su estrecha relación como parte del 
ciclo hidrológico del área. El resto de cuerpos hídricos con datos isotópicos 
disponibles sufre procesos evaporación, lo cual fracciona ambos isótopos 
hacia valores más pesados (Horita, 2005). Esto se observa en las 
muestras Ejido Morelos, Iglesia Nacimiento y Casa Blanca (Punto 2 en la 
Figura 8), así como en las muestras Noria, Ojito y Divi (punto 3 en la 
Figura 8). A modo de resumen se puede evidenciar que los procesos de 
lixiviación y disolución de rocas minerales, la evaporación superficial y el 
transporte de solutos incrementan los contenidos en SOa”-, Cl”, Na* y otros 
cationes, al igual que en los SDT y la conductividad eléctrica en todas las 


muestras en estudio. 


Los datos aquí presentados muestran la complejidad hidroquímica 
de los recursos hídricos en subcuencas semiáridas, lo cual denota que los 
estudios geohidrológicos requieren de un análisis profundo, tanto químico 
como isotópico, para esclarecer el origen de las aguas y su quimismo. El 
estudio también resalta la sinergia que existe entre la hidroquímica y la 
isotopía en la identificación de fenómenos que afectan los cuerpos hídricos 
y, por ende, sirve como caso de estudio para futuros trabajos de este tipo 


en el estado de Coahuila o fuera de éste. 
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Conclusiones 


A partir de los datos y argumentos presentados en este trabajo se 


concluye lo siguiente: 


1. En la zona de estudio se identificaron diferentes tipos de aguas que 
responden a variados procesos geohidrológicos y ambientales, por 


ejemplo: 


a)Se identifican aguas de tipo bicarbonatada-cálcica para los sitios 
129, 143, 13, 137, 135, Ojito, Parquecito y Kikapoo, cuya 
composición refleja la disolución de carbonatos de las formaciones 
geológicas presentes en la sierra de Santa Rosa y sierras 


circundantes. 


b) Las aguas de tipo sulfatada-cálcica y/o magnésica se reportan para 
la noria 17, cuyo enriquecimiento en sulfato sugiere la ocurrencia 
de procesos de intrusión salina o disolución de rocas evaporíticas, 
como el yeso. Dicho mineral podría estar presente en la formación 


La Mula, que aflora en la sierra de Santa Rosa. 


c) Se identifican procesos de salinización en los pozos LF, 124 y CB, 
y las norias IN y EM. El incremento de las sales parece estar 
relacionado con la evaporación superficial y posterior infiltración 
de sales al acuífero. Sin embargo, no se descarta que las 
formaciones de origen delataico, estuarino o tipo sabkha presentes 
en el área pudieran contener minerales solubles como la halita con 


incidencia en el acuífero. 
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2. El diagrama de Mifflin mostró la existencia de diferentes flujos, desde 
local hasta regional. Los flujos locales corresponden con las aguas de 
tipo bicarbonatadas cálcicas; los flujos intermedios son aquellas 
aguas clasificadas como sulfatadas cálcicas, cloruradas sódicas y 
bicarbonatadas sódicas; mientras que las aguas del pozo 124 se 
consideran dentro del flujo regional, pues es la muestra más 
enriquecida en sodio y corresponde a una zona de depositación de 
solutos, donde la composición química tiende a ir aumentando su 
mineralización hasta irse saturando en los diferentes ¡ones a lo largo 
de la trayectoria recorrida por el subsuelo. 

3. Desde el punto de vista isotópico se propone un origen meteórico 
para la mayoría de los cuerpos hídricos en estudio, cuya zona de 
recarga proviene de la sierra de Santa Rosa y las sierras que bordean 
la subcuenta. El resto de las aguas sufren evaporación superficial y/o 
disolución de sales, generando fraccionamiento isotópico y firmas 
isotópicas más pesadas. 

4. Las características geológicas, los ambientes de sedimentación y 
procesos ambientales que ocurren en subcuencas de climas 
semiáridos son determinantes en el variable quimismo observado en 


las aguas y en el deterioro de la calidad del recurso. 
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